Pairing and polarization in systems with retarded interactions 
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In a system where a boson (e.g, a phonon) of finite frequency ljq is coupled to electrons, two 
phenomena occur as the coupling is increased: electron pairing and polarization of the boson field. 
Within a path integral formalism and a Dynamical Mean-Field approach, we introduce ad hoc distri- 
bution function which allow us to pinpoint the two effects. When u)q is smaller than the bandwidth 
D, pairing and polarization occur for fairly similar couplings for all considered temperatures. When 
ujo > D, the two phenomena tend to coincide only for T 3> uio, but are no longer tied for low 
temperatures so that a state of paired particles without finite polarization is stabilized. 

PACS numbers: 71.38.-k, 71.30.+h, 71.10.Fd 



Boson-mediated attractive interactions are known to 
drive instabilities toward broken symmetry phases like 
charge ordering and superconductivity in electronic sys- 
tems. The crossover from weak to strong coupling has 
been extensively studied in both superconducting |l| and 
charge ordered phases 0, 0] , uncovering a general weak- 
to-strong coupling scenario in which at weak coupling 
pairing and condensation occur at the same temperature 
scale while at strong coupling pairing occurs at a higher 
temperature than condensation, leading to an intermedi- 
ate phase of incoherent pairs. 

The same evolution from weak to strong coupling in 
the normal phase has been studied mainly in the limit 
of instantaneous attraction (attractive Hubbard model), 
where the above pairs have been proposed to describe the 
pseudogap regime of the cuprates. This issue, relevant 
also for other field of applications such as the negative 
U localizing centers 0, is in its nature quite hard to 
attack. In fact, the intermediate regime, in which the 
system changes its nature from a Fermi liquid to a paired 
state, intrinsically requires non-perturbative approaches. 

The Dynamical Mean Field Theory (DMFT) emerged 
in the last decade as one of the most successful method 
to deal with non-perturbative regimes. In analogy with 
classical mean-field, DMFT obtains an exact solution for 
local quantum dynamics at the expense of spatial fluc- 
tuations which are frozen, and becomes exact for large 
coordination. Q. DMFT has been successfully applied 
to the study of the Mott-Hubbard transition in the repul- 
sive Hubbard model, in which by increasing the Coulomb 
repulsion, a metal becomes a Mott insulator formed by 
a collection of localized spins. Since the half-filled re- 
pulsive and attractive Hubbard model can be mapped 
one onto the other by a unitary particle-hole transforma- 
tion on one spin species, the Mott-Hubbard transition 
maps onto an analogous metal-insulator transition of the 
attractive model, where the insulating phase is made of 
local pairs and empty sites. This means that much of the 



knowledge we have on the Mott transition can be used 
for the study of the pairing transition of the attractive 
model 0, El • More precisely, we know that the transi- 
tion is a first order line in the interaction-temperature 
diagram, ending in a finite temperature critical point 0] . 
For temperatures larger than the critical temperature, a 
crossover between two phases with different properties 
survives. 

This paper is mainly devoted to discuss how such a sce- 
nario in enriched when the attraction (or equivalently the 
repulsion) has a finite intrinsic energy scale, given by the 
bare frequency of the boson which mediates the attrac- 
tion. The pairing process, which becomes a bipolaronic 
transition when the mediator boson is a phonon, survives 
the inclusion of boson dynamics, and it has been studied 
within DMFT by various authors 0, 0] • The new energy 
scale given by the intrinsic frequency for the boson dy- 
namics has to be compared with both the temperature 
and the electron half-bandwidth D. The most interesting 
region, in which the use of DMFT is crucial, is the inter- 
mediate region in which none of the energy scale and the 
temperature are negligible. The ratio 7 = luq/D between 
the boson energy and the electron hopping energy natu- 
rally defines an adiabatic regime in which 7 is small and 
the electronic degrees of freedom are much faster than 
bosonic ones. In this regime a gradual increase of the 
coupling determines a crossover from quasi-free electrons 
to almost localized "polarons" , quasiparticles in which 
the electrons are strongly bound to the bosonic degrees 
of freedom thereby determining a polarization of the bo- 
son field which couples to the electron density. Such a 
crossover can be unambiguously defined by looking at 
the polarization properties of the boson field 
the opposite antiadiabatic limit, in which the boson dy- 
namics is faster than the electronic one, the interaction 
becomes instantaneous, and we recover either the attrac- 
tive or repulsive Hubbard models. 

The main result of this work is to clearly disentangle 
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the polarization and the pairing crossovers in the nor- 
mal phase and to define unambiguously also the pairing 
crossover. We find that, the larger is gamma, the more 
the two processes occur at different coupling, opening an 
interesting regime in between. 

Although we do not need to specify the physical origin 
of the bosonic mediator, we consider for the sake of defi- 
niteness the Holstein model which describes tight-binding 
electrons coupled to dispersionless Einstein phonons. 
Within DMFT the lattice Holstein model is mapped onto 
a local impurity problem whose action reads [3: 

S„ = - f dr / *X4(«'('-*M (i) 



S el . ph = -VU / drx(T) (n(r) - 1) (3) 
Jo 

where n(r) ~ cI( t ) c t( t )- We take k = mio\ = 1, 
m being the mass of the oscillator. In Eq. Gq 1 ^) 
is the bath propagator that has to be evaluated self- 
consistently 4] . The self-consistency condition which en- 
forces the DMFT contains the informations about the 
original lattice, and for the infinite coordination Bcthc 
lattice we consider it takes the simple form Gq 1 ^) — 
—d T + (D 2 /4)G(r) where G(r) is the local Green's func- 
tion and D is the half-bandwidth. For our model with 
dispersionless Einstein modes, the phononic degrees of 
freedom are not subject to an analogous self-consistency. 
Since the action is manifestly particle-hole symmetric, 
it describes a half-filled system with density n = 1. 
In the antiadiabatic limit (loq — oo, U finite) the ki- 
netic term of the boson field vanishes, the interaction 
becomes therefore instantaneous and the action coin- 
cides with that of an attractive Hubbard model after 
an Hubbard-Stratonovich (HS) decoupling of the quar- 
tic interaction term. We notice incidentally that un- 
der the transformation which maps the attractive Hub- 
bard model on the repulsive one, © becomes S e i^ p h — 
— \/U Jq (1tx(t) (ti^(t) — jt-x(t)) which, once the spin in- 
dex is interpreted as an orbital one, has been studied 
in Ref. El to describe the Jahn- Teller (JT) coupling in 
manganites. 

An efficient method to solve Illl2l.'il) at finite temper- 
ature is the Quantum Monte Carlo method (QMC) in 
the Blankenbecler-Scalapino-Sugar (BSS) approach |12j. 
This method works well in the range of temperatures we 
are interested in (i.e. not too low temperatures) and nat- 
urally yields electronic and bosonic correlation functions, 
as well as the probability distributions associated to the 
boson fields. The method is not affected by the negative 



sign-problem, and its main limitation comes in the adi- 
abatic regime (7 < 1) where the boson becomes heavy 
making more difficult to sample correctly the available 
phase space. In the BSS scheme the fermions are in- 
tegrated out , and the bosons coordinates x(t) are dis- 
cretized into L imaginary-time slices of width At = f3/L 
and then sampled by the Monte Carlo simulation. L has 
to be chosen large enough to reduce as much as possible 
Ar, which controls the the Trotter discretization error. 
To keep Ar less than 1/8 we used 32 slices except for 
the lowest temperature ((3 = 8) for which we have used 
L = 64, As customary, the statistical error is reduced by 
dividing measurements in bins. The self-consistency is 
achieved in a few DMFT iterations, typically from 5 to 
10, depending on the values of the parameters^. 

Now we introduce the key quantities we employ for 
the characterization of pairing and polarization and their 
temperature dependence. Namely, we compute the prob- 
ability distributions of the endpoint 

P(X) = {S(X-x(0))), (4) 

and of the center of mass X c ( "centroid" ) of the boson 
path in imaginary time 




where the averages are evaluated over the action given 
by Eqs. (HEJ. 

The first quantity has a rather straightforward mean- 
ing as a probability distribution of the lattice displace- 
ment (if the boson is thought as a phonon). If P(X) 
has a single maximum, corresponding to a uniform dis- 
placement of the field, the system is not polarized. A 
lattice polarization reflects in the presence of two max- 
ima in P(X), corresponding to opposite polarization of 
occupied and unoccupied sites (bimodal behavior) [ldj . 
In this way a qualitative difference is identified between 
the polarized and unpolarized regimes, which allows for 
an unambiguous way to draw a crossover line, as opposed 
to estimates based on smoothly varying functions as av- 
erage lattice fluctuations or electron kinetic energy. 

The meaning of the centroid variable Xc has been dis- 
cussed in ^3 for a single particle in a binding potential. 
Here the variable X represents the position of the parti- 
cle, and X c is the classical position of the particle [l3| . 
For an heavy particle, the classical limit holds, so that (0} 
and (JSJ coincide |l3| . As the particle is made lighter, the 
wave function becomes broader, increasing the variance 
of P(X) while P(X C ) turns out to be essentially deter- 
mined by the binding range of the potential. Here, we 
use P(Xc) for the many-body problem, and propose that 
pairing, i.e., the formation of local pairs, can be associ- 
ated with a development of a multimodal behavior in our 
centroid distribution P(X C ). This estimator of the pair- 
ing crossover has the same advantage of the previous one, 
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and it finally allows us to draw, beside the polarization 
line U po i, an equally unambiguous pairing line U pa i r (T) 
for any value of the boson frequency. A crossover line 
can also be associated with maxima of susceptibilities re- 
lated to the relevant quantity, e.g., the electron-phonon 
correlation function in the case of polaron crossover [l4| . 
However, this choice is not obvious in the case of pairing, 
where charge, superconductivity or Pauli suscepti- 
bilities can play this role. We defer a comparison with 
these methods in the concluding remarks. 

The ability of our estimator to determine the pair- 
ing crossover can be understood by considering directly 
the interaction term @ in the action. In the adia- 
batic limit (7 — > 0) the kinetic term forces the boson 
path to be r-independent. The boson field becomes 
classical and the interaction term can be written as 
—s/UX Jq dr(n(r) — 1), where X indicates the constant 
value assumed by x(t) along the whole imaginary time 
path. In this limit the centroid coordinate X c is equal 
to X and the two distributions P(X) and P(X C ) obvi- 
ously coincide. Thus the centroid distribution becomes 
bimodal when the system is polarized, which is exactly 
what one expects in the static limit, since a static field 
can induce pairing only with a finite polarization. 

On the other hand, in the opposite atomic (D — > 0, 
7—5-00) limit the electron density becomes a constant 
of motion. Therefore Eq. © takes the transparent form 
— VU(n — 1) Jq <1tx(t) where the electron density is di- 
rectly coupled to the centroid X c . The average appearing 
in Eq. (JSJ) is readily carried out, giving 



P{X C ) (x exp 



,xl 



~ -q log(2 cosh(0VUX c + 1))) 

(6) 

which becomes bimodal for T < U/2. This is exactly 
the scale where double occupancies start to proliferate 
in the atomic limit. Therefore the bimodality of P{X C ) 
correctly signals the onset of pairing also in the antiadia- 
batic regime. In the same limit, it can be proved that the 
endpoint distribution P(X) has a variance which scales 
with 1/VAr and as a consequence no definite polariza- 
tion may occur. We finally notice that the D — > limit 
of adiabatic P(X) 10] coincides with P{X C ) of Eq. JBJ. 
Since for ui — the distributions of X and X c coincide, 
we conclude that in the atomic limit P(X C ) is the same 
for ujq = and ujq = 00. This suggests that the pair- 
ing crossover may depend on loq more weakly than the 
polarization one. 

To analyze the evolution of P{X) and P(X C ) at finite 
D and ujo we use BSS-QMC. The numerically exact re- 
sults, shown in FigQ] clearly show that P(X) and P (X c ) 
tend to coincide in the relatively adiabatic case 7 = 0.1, 
as expected from the previous arguments about the adi- 
abatic limit. The two quantities are clearly different for 
7 = 1 and 8. Actually for temperatures smaller than u>o, 
the polarization crossover U po i moves to larger values as 



7 is increased, while the line (U pa i r ) where P(X C ) be- 
comes bimodal is only slightly shifted to larger couplings 
with increasing 7. This is strongly reminiscent of the be- 
havior of the metal-insulator transition in the Holstein 
model at T = 0, whose critical coupling is slowly increas- 
ing with 7 and then saturates to the asymptotic 7 = 00 
value 0- On the other hand, the polarization crossover is 
almost directly proportional to 7pJ. Both at zero and at 
finite temperature the line where the centroid becomes 
bimodal does not coincide with the metal-insulator line, 
but it can be considered a precursor which depends in a 
very similar way on 7 and on U / D. 
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FIG. 1: (color online) Behavior of U pa ir{T) (in blue and with 
the letter c) and U po i(T) (in red and with the letter x) at 
7 = 0.1, 1 and 8, respectively indicated by solid and open 
triangles, solid and open circles and solid and open squares. 
The green dashed line represents the 7 = result for both 
Upair (T) and U P oi(T). The dashed arrows indicates the zero- 
temperature results for the polaron crossover 0. 



Our DMFT results can also be compared with the 
semi-analytical results for 7 = [lCj, which is repre- 
sented by the green dashed line in Fig^ It can be seen 
that the 7 = 0.1 case is in very good agreement with the 
adiabatic result, and also the 7 = 1 and 8 cases, at high 
temperature, fall on the same curve. It must be observed 
also that, in the full DMFT calculation, the centroid 
distribution depends weakly on 7, as suggested by the 
atomic limit. Interestingly, the adiabatic result displays 
a re-entrance at low temperatures. Although the QMC 
simulations do not reach sufficiently low temperatures, 
we find that the re-entrance is present also for 7 differ- 
ent from zero, as indicated by the arrows in Fig^ which 
mark T = results for the polarization crossover for 
the Holstein modelQ- This phenomenon has been also 
reported for the Mott transition in the repulsive model 
|l6j, and associated to spin entropy of the insulator. The 
same physics holds here, where the entropy is associated 
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to a pseudo-spin. 

A strong indication that the centroid distribution 
marks the pairing crossover is provided by the analy- 
sis of the double occupancies (n^n±). This quantity in- 
creases from the non-interacting value 1 /4 to the asymp- 
totic value of 1/2 as the coupling is increased. We com- 
puted (n^rii) for different temperatures in the antiadia- 
batic regime 7 = 8, chosen in order to have clearly differ- 
ent values for U po i and U pa i r - We notice that our results 
are correctly close to the attractive Hubbard model. 
In Fig[21 Upoi and U pa i r are marked, respectively, by a 
black solid square and a blue solid circle, for each temper- 
ature. {n-\n\) at U po i is rapidly increases with decreasing 
temperature, while (n^n;) calculated at U pa i r is almost 
constant with temperature, or even slightly decreasing 
at low T. So, when the centroid distribution becomes 
multimodal, the double occupancies are approximately 
70% of the saturation value of 1/2. Similar results are 
found also for 7=1 and 0.1, where it is however harder 
to distinguish U pa i r from U po i- This establishes the ef- 
ficiency of our centroid distribution in pinpointing the 
pairing crossover for the whole range of parameters. 
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FIG. 2: (color online) Double occupancies for different tem- 
peratures at 7 = 8. The blue solid circle (black solid square) 
denotes the coupling where P(X C ) (P(X)) develops the bi- 
modal character. 

In conclusion for a half-filled system subject to re- 
tarded attractive (repulsive) interaction we observe three 
qualitatively different behaviors i) a normal state of dclo- 
calized particles, ii) a pair region of bounded incoherent 
pairs with no associated polarization, in) a bipolaron re- 
gion in which pairing is associated with a definite polar- 
ization. These three regions are delimited by two clear- 
cut lines where the probability distributions for the cen- 
troid and for the endpoint of the bosonic path, develop 
a bimodal shape. While the border between the nor- 
mal and the local pair phase is weakly dependent on the 



adiabatic ratio, the one between the pair and the bipo- 
laronic phase, shifts to higher value of the coupling with 
increasing 7, so that the pair phase gets larger, even- 
tually covering the whole right hand part of the phase 
diagram in the attractive Hubbard limit (7 = 00). The 
Upair line is very similar to the T* line found in Ref.|l5| 
by looking at the maximum of the spin susceptibility. For 
temperatures of the order of u>o the two lines merge be- 
cause phonons at such temperatures behaves classically. 
Finally, due to the electron-hole symmetry the retarded 
attractive interaction is mapped onto a JT model leading 
to an Hubbard- like repulsion in the antiadiabatic limit. 
In this case the crossover at U pa i r is associated to the 
formation of local orbital ordering for JT interactions. 
The crossover at U po i survives only in the retarded JT 
case and is shifted toward higher value of the coupling as 
the phonon frequency increase leaving a region in which 
orbital ordering is not associated with JT polarons. 

We acknowledge useful discussions with C. Castellani 
and A. Toschi. This work was supported by MIUR-Cofin 
2003 matching funds programs. 
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